use "[1] Flooding map data.dta"

grmap mean_to_city_ratio using china_city_long4490_coord.dta, id(ID) fcolor("252 255 255" "224 224 224" "176 176 176" "128 128 128" "80 80 80" "32 32 32" "0 0 0" ) ocolor("black" ...) clmethod(custom) clbreaks(0 1.17 2.34 3.51 4.68 5.85 7.02 1000) graphr(margin(medium)) osize(*0.01 ...) legend(size(*1.0) pos(7) order(2 "0 to 1.17%" 3 "1.17% to 2.34%" 4 "2.34% to 3.51%" 5 "3.51% to 4.68%" 6 "4.68% to 5.85%" 7 "5.85% to 7.02%" 8 "7.02% and above") ti(Sum, size(*0.5) pos(11) color(black)) color(black)) line(data(china_city_long4490_line_coord.dta) size(*0.5 ...) color(black) select(keep if _ID <= 5)) label(data(china_city_long4490_label) x(x_centroid) y(y_centroid) l(label) color(black) size(*0.8)) xsize(10) ysize(10) aspectratio(0.9) freestyle ysc(off) xsc(off)
